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Abstract 

Calogero-Moser models can be generalised for all of the finite reflection groups. 
These include models based on non-crystallographic root systems, that is the root sys- 
tems of the finite reflection groups, H^, H4, and the dihedral group him), besides the 
well-known ones based on crystallographic root systems, namely those associated with 
Lie algebras. Universal Lax pair operators for all of the generalised Calogero-Moser 
models and for any choices of the potentials are constructed as linear combinations of 
the reflection operators. The consistency conditions are reduced to functional equa- 
tions for the coefficient functions of the reflection operators in the Lax pair. There are 
only four types of such functional equations corresponding to the two-dimensional sub- 
root systems, A2, B2, G2, and l2{m). The root type and the minimal type Lax pairs, 
derived in our previous papers, are given as the simplest representations. The spec- 
tral parameter dependence plays an important role in the Lax pair operators, which 
bear a strong resemblance to the Dunkl operators, a powerful tool for solving quantum 
Calogero-Moser models. 



1 Introduction 



Generalised Calogero-Moser models are integrable many-particle dynamical systems based 
on finite reflection groups. Finite reflection groups include the dihedral groups /2(^) and H-^ 
and Hi together with the Weyl groups of the root systems associated with Lie algebras, called 
crystallographic root systems. Integrability of classical Calogero-Moser models based on the 
crystallographic root systems |Jl|, |^ is shown in terms of Lax pairs. The root and the minimal 
type Lax pairs derived in our previous papers 0] provide a universal framework for these 
Calogero-Moser models, including those based on exceptional root systems and the twisted 
models. On the other hand, a theory of classical integrability for the models based on non- 
crystallographic root systems has been virtually non-existent. This is in sharp contrast with 
the quantum counterpart. Dunkl operators, which are useful for solving quantum Calogero- 
Moser models, were first explicitly constructed for the models based on the dihedral groups 

i- 

In this paper we present a Lax pair in an operator form for generalised Calogero-Moser 
models, which applies universally to the models based on non-crystallographic root systems 
as well as those based on crystallographic ones. This Lax pair, as expected, bears a strong 
resemblance to the Dunkl operators ^ and the reflection operators play a central role. 
The spectral parameter dependence is also essential, in contradistinction to the auxiliary 
role in the conventional formalism 0. When suitable representation spaces are chosen, the 
universal Lax pair reproduces the root type and the minimal type Lax pairs for the models 
based on the crystallographic root systems. This provides another clue that the quantum 
and classical integrability of the generalised Calogero-Moser models are closely connected. 
We hope that this Lax pair operator formalism will cross-fertilise the fruitful subject of the 
quantum and classical Calogero-Moser systems and related ones such as Toda systems. 

For the general background and the motivations of this paper and the physical appli- 
cations of the Calogero-Moser models with various potentials to lower-dimensional physics, 
ranging from solid state to particle physics and supersymmetric gauge theories, we refer to 
our previous papers and references therein. 

This paper is organised as follows. In section |^ we summarise the basic concepts of the 
finite reflection groups in order to set the stage and to introduce appropriate notation. In 
section ^ the generalised Calogero-Moser models are defined with various choices of root 
systems and potentials. Section ^ is the main body of the paper. The Lax operators L and 
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M are defined as a linear combination of the reflection operators Sp for all the roots p. The 
coefficient functions depend on the dynamical variables and on the spectral parameter in a 
rather symmetrical way. Consistency of the Lax pair can be proved quite easily. The Lax 
equation is at most quadratic in the reflection operators, SpSo-, which becomes an identity 
operator for p = a, = 1, and a two-dimensional rotation operator for p ^ a . The linear 
(s^) and the constant parts give the canonical equations of motion of the generalised 
Calogero-Moser models. The quadratic part SpS^^ P 7^ cr, imposes the consistency conditions, 
which are decomposed into those corresponding to two-dimensional sub-root systems con- 
taining p and a. As shown in Table 1, there are only four types of two-dimensional root 
systems, A2, B2, G2, and /2(5) for all of the Coxeter groups except for the dihedral group 
/2(m) which can have many for some values of m. Thus the functions appearing in the Lax 
pair operator (except for those for l2{m)) need to satisfy at most two functional equations. 
The solutions are derived in the Appendix. In section |^ various possible representations of 
the Lax pair operator are discussed. The minimal type and the root type Lax pairs are 
derived as two simplest examples in section ^. Various sum rules utilised in previous papers 
are derived as restrictions of the general functional equations derived in section ^. Some 
comments and discussion are given in section |^. Details of the derivation of the solutions are 
relegated to the Appendix. Two types of solutions, the untwisted and twisted solutions, are 
derived. The consistency conditions of all of the untwisted solutions are attributed to one 
functional identity, ( [A.12| ), of the Weierstrass cr function. Twisted solutions are obtained as 



proper linear combinations of the untwisted ones. 

2 Root systems and finite refiection groups 

We now review some facts about root systems and their reflection groups in order to introduce 
notation |^. We consider only reflections in Euclidean space. A root system A of rank r is 
a set of vectors in R*" which is invariant under reflections in the hyperplane perpendicular 
to each vector in A. In other words, 

s«(/3)gA, Va,/3GA, (2.1) 

in which 

s^{(3) = (3 -2{a- (3/\a\'^)a. (2.2) 
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Dual roots are defined by = 2q;/|q;P, in terms of which 

SaiP) = /? - (a^- P)a. (2.3) 

The orbit of /9 G A is the set of root vectors resuhing from the action of the reflections on it 
{sa{(3),C( G A}. The set of positive roots A+ may be defined in terms of a vector V G R*", 
with • q; 7^ 0, Va G A, as those roots a G A such that a ■ V > 0. A unique set of r 
simple roots 11 is defined such that they span the root space and the coefficients {aj} in 
P — ttjaj for (3 G A+, {aj & U, j — 1, . . . ,r} are all positive. 

The set of reflections {sq, q; G A} generates a group, known as a Coxeter group. It is 
generated by products of Sa with a G H subject to the relations 

(saS/^r^"'^^ = 1, a,(3eU. (2.4) 

The set of positive integers m{a,pi) uniquely specify the Coxeter group with m{a,a.) — 
1, Va G n. For example, for two-dimensional crystallographic root systems ^2, B2, and G2, 
the integer m{a,f3) is 3, 4, and 6, respectively. Thus Sq,S/3 is a two-dimensional rotation by 
an angle ±27r/3, ±7r/2, and ±7r/3, respectively. This fact will be used in later sections. We 
consider here only those Coxeter groups with a finite number of roots in Euclidean space, 
called the finite refiection groups. 

The root systems for finite refiection groups may be divided into two types: crystallo- 
graphic and non-crystallographic root systems. Crystallographic root systems satisfy the 
additional condition 

a'^-f3eZ, Vq;,/3gA. (2.5) 

These root systems are associated with simple Lie algebras: {Ar, r > 1}, {Br, r > 2}, 
{Cr, r > 2}, {Dr, r > 4}, Eg, E7, Eg, F4, and G2 and {BCr, r > 2}. The Coxeter groups 
for these root systems are called Weyl groups. The remaining non-crystallographic root 
systems are H^, H4, and the dihedral group of order 2m, {l2{m), m > 4}. 

Weyl chambers are defined as the open subsets of R'' that result from removing the 
refiection hyperplanes Ha, a G A, if^ = G R^, a^- q — 0}. The action of the refiection 
group on the Weyl chambers is transitive and free, i. e. any Weyl chamber may be obtained 
from another by the action of an element of the refiection group and this element is unique. 
The principal Weyl chamber is defined as the one whose points q satisfy q ■ a > 0, Va G A+. 
For crystallographic root systems this implies that all points in the principal Weyl chamber 
have positive Dynkin labels. 
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Definitions and properties of tlie crystallographic root systems may be found in many 
references, see for example 0. The concept of weights may be defined for these root systems. 
For crystallographic root systems, any positive root is a sum of simple roots in 11 with 
positive integer coefficients. The set of weights S, which lie on a lattice, is the set of 
vectors such that if A G S then a^- A is an integer for any a G A. Fundamental weights 
are vectors which form a dual basis to the corresponding dual simple roots a^, i.e. 
aj ■ A*^'^') = j, k = 1, . . . ,r. Any weight \E' G S may be expressed as a sum of fundamental 
weights with integer coefficients, \E' = I]j=i '^j^^"'''- The coefficients {dj, j = l,...,r} are 
called the Dynkin labels of 

We now briefly describe the non-crystallographic root systems. The dihedral group of 
order 2m, l2{fn), is the group of orthogonal transformations that preserve a regular m-sided 
polygon in two dimensions. It consists of m rotations (through multiples of 27r/m) and m 
reflections. The angle between adjacent roots is vr/m and a possible basis for the roots, if all 
are chosen to have the same length jajp = 1, is 

aj = (cos(j7r/m), sin(j7r/m)) , j = 1, . . . , 2m. (2.6) 

For odd m all of the roots are in the same orbit of the reflection group but for even m there 
are two orbits, one consisting of the aj with odd j and the other with even j . This then 
allows two different coupling constants and potential functions for the l2{m) Calogero-Moser 
model for even m. 

The reflection group of type is the symmetry group of a regular 120-side solid, with 
dodecahedral faces, in R^. It is a group of order 14400. The group of type is a subgroup 
of and is the symmetry group of the icosahedron (with 20 faces) in R^. It is a group of 
order 120. Deflne 

TT l + y/S 2n -1 + V5 

a = cos— = : , o = cos— = . 2.7) 

5 4 5 4 

Then a choice of simple roots for is the following: 

ai = (^a,~,b,0^ , a2 = (^-a,^,b,0^ , (2.8) 
as = Q,6, -a,0^ , = (^-^, -a,0,b 

The full set of roots of if4 in this basis may be obtained from (1,0,0,0), (|,|,|,|), and 
(a, |, 6, 0) by even permutations and arbitrary sign changes of coordinates. These 120 roots 



form a single orbit. A subset of (^78|) , {ai, 02,0^3} is a choice of simple roots for the H3 
root system. In this basis, the full set of roots for results from even permutations and 
arbitrary sign changes of (1, 0, 0) and (a, ^,b). These 30 roots also form a single orbit. 

3 Generalised Calogero-Moser Models 

A generalised Calogero-Moser model is a Hamiltonian system associated with a root system 
A of rank r. Quantum versions of these models are also integrable, at least for certain choices 
of A and potential function [Q. The dynamical variables are the coordinates {q^} and their 
canonically conjugate momenta {pj}, with the Poisson brackets 

{(l\Pk} = Si, {q^,qk} = {Pj,Pk} = 0, j,k = l,...,r. (3.1) 

These will be denoted by vectors in R'' 

g= p = (pi, . . . (3.2) 

The Hamiltonian for the Calogero-Moser model is 

n = lp'+j:riV\^\{a-q), (3.3) 

in which the real coupling constants g\a\ and potential functions V\a\ are defined on orbits of 
the corresponding finite reflection group, i.e. they are identical for roots in the same orbit. 
This then ensures that the Hamiltonian is invariant under reflections of the phase space 
variables about a hyperplane perpendicular to any root 

Q Sa{q), P Sa{p), Va e A (3.4) 

with Sa defined by (|2.3|) . 

The Lax pair operator that we will construct in later sections will apply for the following 
potentials {g\a\ = g for all roots in simply laced models and g\a\ = gi for long roots and 
g\a\ = Qs for short roots in non-simply laced models): 

1. Untwisted elliptic potential. This applies to all of the crystallographic root systems 
and the potential function is 

V\a\{a ■ q) = p{a ■ q\{2uJi,2uj^}), for all roots, (3.5) 
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in which p is the Weierstrass p function with a pair of primitive periods {2uji,2uj^}. 
Throughout this paper we adopt the convention that the Weierstrass p, C,, and a 
functions have the above standard periods, unless otherwise stated. 

2. Twisted elliptic potential. This applies to all of the non-simply laced root systems. 
Except for the G2 model, the potential functions are 

{p{a ■ g|{2a;i, 2a;3}), for long roots, 
(3.6) 
p{a ■ q\{jJi, 2a;3}), for short roots. 

That is, the potential for short roots has one half of the standard period in one direction, 
which we choose to be uj^. For the G2 model, 

' p{a ■ q\{2ijji, 2a;3}), for long roots. 



V\a\{a ■q) = < 



(3.7) 

^ p{a ■ q\{^, 2a;3}), for short roots. 

In this case the potential for short roots has one third of the standard period in one 
direction, which we choose to be (Vi. The cases of BCr system and the extended B^., 
Cr, and BCr systems will be discussed separately in later sections. 

3. Trigonometric and hyperbolic potentials. This applies to all crystallographic systems 
and the potential functions are 

{a^ /siv? a{a ■ q), 
for all roots, (3.8) 
a^/sinh^ a{a ■ q), 

in which a is an arbitrary real constant. 

4. Rational potential. This applies to all of the generahsed Calogero-Moser models includ- 
ing those based on the dihedral group hiiTi), H^, and H4 and the potential function 
is 

V|Q,|(a • q) = ^, for all roots. (3.9) 

These models are also integrable if a confining harmonic potential 

(3.10) 

is added to the Hamiltonian. 
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Some remarks are in order, (i) For all of the root systems and for any choice of the 



1 



potential, the generalised Calogero-Moser model has a hard repulsive potential ~ Jl^^ 
near the reflection hyperplane = {q & R*^, a^- q = 0}. This repulsion potential is 
classically insurmountable. Thus the motion is always confined within one Weyl chamber. 
In other words, the spatial ordering of the particles is unchanged during the time evolution. 
This simplifying feature is a basic cornerstone of the solvability. The coupling constants 
(with a scale l/lap) are measures of the strength of the repulsive potentials, (ii) The 
trigonometric, hyperbolic, and the rational potentials are obtained from the elliptic potential 
as one or both periods tend to infinity. The Lax pairs for these degenerate potentials can 
be obtained from the one for the elliptic potential by taking the corresponding limit. Thus 
we do not write down the Lax pairs for degenerate potential cases except for the models 
based on the non-crystallographic root systems, for which only the rational potentials are 
integrable. (iii) For all of the Lax pairs based on any root systems and any choice of the 
potential, except for the rational potential with the harmonic confining potential (|3.10|) , one 
can introduce an additional complex parameter ^ {spectral parameter), which appears in 
the equation for the spectral curve |10|- (iv) Independent conserved quantities Tr{L^) 



to be obtained from a Lax equation L = [L, M] occur at such k = 1 + exponent of the 
corresponding crystallographic root systems. For the non-crystallographic root systems, 
they arise at k = 2,m for the dihedral group l2{ni), k = 2,6, 10 for and k = 2, 12, 20, 30 
for H4. These are the degrees at which Coxeter invariant polynomials exist [Q. 

4 Lax pair and functional equations 

Here we construct a Lax pair for the generalised Calogero-Moser models in an operator 
form acting on an as yet unspecified vector space and derive the necessary and sufficient 
conditions for the consistency of the Lax equations. The Lax pair ( [4.5D contains operators 
as well as functions x\p\{u, w), y\p\{u, w), which are related to the chosen potential (|3.5|) - (|3.9|) . 
The consistency of the Lax pair requires that x\p\{u,w) satisfy certain functional equations 
(|4.27| ), ( |A.7|) , which are closely related with those required for the commutativity of Dunkl 



operators |, §. Verification that the solutions j ^Ml , (H), ( CT) , (HID, (H), (EH) 



satisfy the functional equations will be presented in the Appendix. 

The operators appearing in the Lax pair for a generalised Calogero-Moser model asso- 
ciated with a root system A are naturally the reflection operators {sa, a G A} of the root 
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system. They act on a set of R*" vectors T = {jj,^''^ G R*", k = 1, . . .}, which is closed under 
the action of the reflection group. The totality of the vectors in F forms the representation 
space V. A general construction of the representation space, and some explicit cases will 
be presented in the subsequent two sections. Another set of operators {Hj, j = 1, . . . , r} is 
necessary. If Hj acts on a vector n^''^ E T, the j-th component is returned: 

These form the following operator algebra: 

[H,,Hk] = 0, (4.1) 
= a,(a^-^)s,, (4.2) 

{s^rsf.r^'^^f^ = 1. (4.4) 

The first relation ( |4.1| ) implies that the operators {Hj, j = 1, . . . ,r} form an abelian sub- 
algebra and relations ( [4.3| ) and ( [4.4| ) are just those for the finite reflection group associated 
with the root system A. The set of integers m{a,(3) are those appearing in the Coxeter 



relations (|2.4|) which characterise the reflection group. 

Next we describe the Lax pair and the corresponding Hamiltonian for the generalised 
Calogero-Moser model for the root system A. The Lax operators are 

L = p-H + X, X = iY. g\,\ (p^- H) x\,\ (p • g, (p^- H)i)sp, (4.5) 

PGA+ 

M = i 9\p\y\p\{p-(lAp^-H)C)Sp. 

P6A+ 

The function y\p\ is defined by 

d 

y\p\iu,w) = —x\p\{u,w). (4.6) 

The variable u takes care of the dynamical variable dependence and w is for the spectral 
parameter dependence. Furthermore, it is required that x\p\{u,w) is odd: 

x\p\ {-u, -w) = {u, w) (4.7) 

so that L and M are independent of the choice of positive roots A+. This also implies that 
the sums in ( [4.5|) may be extended to a sum over all roots if an additional factor of 1/2 is 
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included in front of the sums since the summands are even under p 
x\p\{u, w) is assumed to have a simple pole at m = with unit residue 

lim ux\o\iu, w) = 1. 



-p. The function 



(4.8) 



This condition is related with the unit strength of the repulsive potential near the reflection 
hyperplane mentioned earlier. 

Before proving the consistency of this Lax pair we calculate the Hamiltonian and find 
the corresponding equations of motion. The Hamiltonian for the theory is defined in terms 
of a representation of the operator L of ( [4.5| ) by 

1 



n 



2Cr 



Tr{L' 



where the constant Cr, which depends on the representation, is defined by 
The resulting Hamiltonian is then 



\P\ 



(4.9) 



(4.10) 



(4.11) 



in which C is independent of the dynamical variables q and p and therefore unimportant for 
the classical theory. The potential functions are given in ( p.5|) -( pl9D . The function x\p\ and 
the potential function V\p\ are related (except for the confining harmonic potential (|3.1CI|) , 
which has to be added separately) simply as 



x\p\{u,w)x\p\{-u,w) = -V\p\{u) + C\p\{w), 



(4.12) 



i.e. the product gives a sum of a function of only u and a function of only w. It is easy 
to show that all of the functions x\p\{u,w) which lead to a consistent Lax equation, (|4.34| ), 
(lOTI ), (|0^ ), ( CT ), ( P^ ), (1^, satisfy this property. 

The equations of motion following from this Hamiltonian are 

m _ 

m 



(4.13) 



d_ 



\P\ 



w 



peA 



(4.14) 
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Because of (|4.12| ) the last expression in ( [4.14| ) is independent of w. 
The Lax equation 

L = [L, M] 

may be divided into three parts as 



(4.15) 



ip-H) = [X,M] 



diagonah 



(4.16) 
(4.17) 

= [X,M]off.d^agonal. (4.18) 

We next prove, in turn, that each of these equations is consistent with the equations of 



d 
di 



motion ([4.13|) . 

The left-hand side of ( tl.l6| ) is 

d 



dt 



^ = ^ 9\p\ {p"- H) y\p\{p ■ q, {p"- H)^) (p ■ q) 



(4.19) 



and the right-hand side is 



[p-H,M] = [p-H,z J2 g\,\y\,\{p-q, (p^- H)0 s 

PGA+ 

= ^ J2 9\p\ y\p\ (p ■ {p^- H)Cj [p-H,s 

PGA+ 

= ^ 9\p\ y\p\{p ■ q, {p^- H)0 (p^- H) (p ■ p) s 

d 



P' 



dt^ 



-X. 



(4.20) 

and the last line follows from the 



The third line follows from the commutation relation 
equation of motion q = p. 

The left-hand side of ( [4.1?| ), after using the equations of motion (|4.13| ) is 



d 9 

-tAp-H) = Y r^ip-^) [y\p\{p ■ q,w)xipi{-p ■ q,w) - xipiip ■ q,w)yipi{-p ■ q,w)] . (4.21) 

pGA \P\ 

The summand x\p\ {u, w)y\p\{—u, w) — x\p\ {—u, w)y\p\ {u, w) is independent of the parameter w. 
The commutator [L, M] reads 

[X,M] = ~[Y: g\p\ {p^-H)x\p\{p-q,{p''-H)i)sp, 9w\yUa ■ q,{a^- H)i)s„\ 

peA+ o-eA+ 

= - Y aipigwilip"^- H)x\p\{p-q,{p'^- H)^)y\^\{(7-q,{sp{(7y- H)^)SpS^ 
p,(TeA+ 

-1/|.|((T . q, (a^- H)0 (s^ipy- H) X|,|(p ■ g, (s.(p)^- H)0 sjp] .(4.22) 
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Since SpSa and s^Sp are rotations (except for p = a, = 1) they do not leave any real vectors 
in the rotation plane invariant. Thus [X, M] is decomposed into the diagonal (p = a) and 
the off-diagonal (p 7^ cr) parts. The diagonal part gives the equation of motion 

[X, MUa,. = E 9Up''- H) [y\p\{p ■ q, (P^- H)i) X|,|(-p ■ g, (p^- H)i) 

PGA+ 

-a;|,|(p ■ g, (p^- y|p|(-p ■ g, (p^- (4-23) 
= Jt^P-H)- (4.24) 

Finally, ([4.1 8|) will lead to the functional equations which must be satisfied by the func- 
tions x\p\{u, w) in the Lax pair. Relabeling the roots in the second summation of (|4.22| ) gives 
the general operator form for the functional equation 

= [X,MU_,,ag, 

= E 9\P\9W\ [(P^- H)x\p\ (p ■ g, (p^- H)i)yUc^ ■ q, (s,(a)^- H)0 

~{sp{ar- H)y\p\{p ■ q, (p^- H)Oxia\{a ■ g, (sp(a)^- H)0]SpS^. (4.25) 

This shows that the consistency condition involving all of the roots is decomposed into 
a sum of two-dimensional ones corresponding to a fixed rotation = SpS^ in each two- 
dimensional plane. Since the coefficient of = SpS„ in this equation must separately 
vanish, this may be decomposed into the functional equations 

= E ^|p|^ki[(p'^-p)a;|p|(p-g, (p'^-p)0z/ki(f^-9, (sp(fT)^-p)0 

-{sp{ay- p)y|p|(p ■ g, {p^- p)0 x\^\{(T ■ g, (sp((T)^- p)0] , (4.26) 

in which p is a generic vector in F. This equation must be satisfied for a fixed rotation 
R^ = SpSrj and all roots appearing in it are in the two-dimensional sub-root system $ = 
{k, k G {^r\span{p, o")} with positive roots = $nA+. The only possible two-dimensional 
root systems $ are Ai x Ai, A2, ^2, and l2{fn)- Table 1 shows the two-dimensional 
sub-root systems appearing in the root systems of finite refiection groups. The Ai x Ai root 
system has been omitted since its corresponding functional equation is trivially satisfied for 
any function and therefore does not give any constraint on the functions x\p\{u, w). It should 
be stressed that the functional equations are determined by the two-dimensional sub-root 
systems only and not by where they are embedded in the entire root system. Thus, as seen 
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Root System 


Sub-root Systems 


Ar, r > 1 


A2 


Br,r>2 




C„ r > 2 




A', ^ > 3 


^2 

-^z 


BCr, r > 2 


^2,^2 


Ee,Ej, Eg 


^2 


^4 


^2,^2 




^2,6-2 


/2(m) 


h{ky 




A2,/2(5) 


^4 


^2,/2(5) 



Table 1: Two-dimensional sub-root systems. Ai x Ai is not included, f: k divides m. 

from Table 1, each function x\p\ must satisfy at most two functional equations, except for 
the models based on the dihedral groups. 



= /X ■ R^/2 — 



The functional equation in ( |4.26| ) may be further simplified to 

d_ 
dq 

in which 

[a-- p) 



(4.27) 



E 



x\p\{p ■ (P^- P)0 ■ g, {sp{ay- p)0 (4-28) 



and -R^/2 is a rotation by an angle iIj/2 with i?^ = SpS^- It has the following action on any 
pair of roots p, a G $+ which satisfy = SpS^ 



D \P\ {<^^- P) I N 

R^I2P = ------ Sp{a), 

\a\ \a^- p\ 

R^/2 o- = — |_,, _| p. 



(4.29) 



Using these relations, the simplified form of the functional equation ( [4.27|) may be proven. 
First substitute the definition of f{q,p) in ( [4.28|) into (|4.27| ) to obtain 

(av. p) - 



J2 9\p\9\a\ 



pJ^a£^+,r=Sp 



W^- p\ 



(/i ■ R^/2 p)y\p\ip -q^P ■ P)0 ■ (sp(c^) ■ p)0 



+ (/i ■ R^/2 (r)xipi{p ■ q, {p"^- /i)0 y\a\{(y ■ q, (sp(a)'^- p)0 



(4.30) 



Since the coupling constants g\a\ are arbitrary and are constants on the orbits they may be 
rescaled as g\a\ 2g\a\/\<y\ for all a G $. Performing this rescaling of the coupling constants 
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and using the action of in ( |4.29| ), the previous equation becomes 

2 



= 51 9\p\g\a\ 



/^)2/|p|(p ■ Q, {p^- P')0 ■ isp{(^y ■ P')0 



. (4.31) 



This is the same as the earher form of the functional equation in ( [4.26|) , after using the 
definition = 2a/\a\'^. 

We next present the functional equations for the two-dimensional root systems, A2, B2, 
G2, and l2{fn), and state their solutions. The proofs that the solutions satisfy the relevant 
functional equations are contained in the Appendix. 

We first consider the functional equation for the A2 root system with simple roots a and 
p. We choose R2n/3 = SaSi3 and the functional equation ( [1.25| ) becomes 

= /i-i?^/3 ^fAMf^) (4.32) 

in which {q, p) is defined by 

fA2{q, P) = x{{a + (3)-q,{a + (3^- /i) x{a ■ g, /x) 
+ ■ g, yu) x((q; + /5) ■ g, a^- /i) 

— x{a ■ q,a'^- fi) x{P ■ q, {a + fi). (4.33) 

The subscripts on the function x is omitted since all roots are in the same orbit and hence 
only one function appears. Also the spectral parameter ^ is absorbed into fi by redefinition: 
fi ^ fi/C,- We look for solutions to this functional equation, as well as the ones for other 
root systems, which are valid for arbitrary vectors q and fi. Therefore these solutions are 
valid in any representations. In certain representations, such as the minimal and the root 
type representations discussed in section H, there is a larger class of solutions. The functional 
equation arising from the rotation R-2n/3 is equivalent to that given above. A simple solution 
to ( [4. 321) satisfying the residue condition ( ^l8|) is 



x(u, w) = '^1^, — ^explbwu], (4.34) 
a{w)a{u) 

in which b is an arbitrary complex parameter. The potential function in the resulting Hamil- 
tonian (|3.3| ), V{u) = p{u), is the Weierstrass elliptic function with the set of primitive 
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periods lujx and 2co'3. The limits as one or both periods diverge yield the potential functions 
a} j sin^(au) and j sinh^(aM) or respectively 0. 

Functional equation (|4.26|) for the root system may be written in terms of the short 
and long simple roots, a and /3, respectively. We choose i?,r/2 = s^sp in ([4.25|) and the 
resulting functional equation is 

d_ 

dq' 

in which 



= /i ■ —fB2 (q, At) (4.35) 



fB2{Q,l^)= - xs{a-q,a'^- n)xL{p-q,{2a + py- n) 

+ xs{{a + (3)-q,{a + (3^- fi) XL{{2a + f3) ■ q, fi) 

+ xs{a ■ q, -{a + P^- /i) XL((2a + P) ■ q, (2a + p^- n) 

+ xs{{a + P) ■ q,a^- ^)xl{P ■ q-,0^- n)- (4.36) 

In this case there are two sets of elliptic function solutions for the long and short root 
functions, Xl and Xs-, respectively. The solutions along with the corresponding potential 
functions are 

xl[u,w) = xs[u,w) = exp[bwu\, (4.37) 

(J[W)(T[U) 



Vl{u) = Vs{u) = p{u) 



for the untwisted solution and 



a{w-u) 
xl{u,w) = exp[bwu\, 
a{w)a{u) 

(t{w/2-u\{uji,2uj3}) ei 

XsiU.W) = — —77 r- —77 r-explOH ) W U\ 

' ' a{w/2\{uji,2u;3})a{u\{ui,2ius}) 2' ^ 

_ xl{u,w/2)xl{u + uji,w/2) 
xl{uji,w/2) 

VLiu) = p{u), Vs{u) = p{u\{uji,2uj3}), 



(4.38) 



for the twisted solution. Here ei = p{uJi) and b is an arbitrary complex constant. 

Next we consider the functional equation ( 4.26| ) for the G2 root system with short and 



long simple roots a and P, respectively. We choose i?^/3 = SaSp in ( |4.26| ) and the functional 
equation is 

= /i-i?^/6 ^/G,(g,/i), (4.39) 
15 



in which 



/g2(?,/^)= - Xs{a-q,a'^-ij)xL{f3-q,{Sa + (3)'^- jj,) 

+ xs{a ■ q, -{2a + p^- fx) ^^((Sa + P) ■ q, (3a + p^- fi) 

+ xs{{2a + p)-q, (2a + p^- /i) ^^((Sa + p) ■ q, -(3a + 2/?)^- fx) 

+ xs{{2a + p)-q,-{a + fx) XiH^a + 2P) ■ q, (3a + 2/?)^- fx) 

+ xsiia + P)-q,{a + fx) XL{{3a + 2/5) ■ q, /i) 

+ a;s((a + /3)-g,a^-/i)a;L(/3-g,/3^-/i). (4.40) 

As before, xs and xl are the functions for short and long roots, respectively. The elliptic 
function solutions to (|4.39|) along with the corresponding potential functions Vs{u) and Vl{u) 
are the untwisted ones in (^4.37|) and the following twisted ones 

a{w-u) 
xl[u,w) = , , , , exp[bwu\, 
a[w)a[u) 

I \ o(w 13 ~ u\{2ijj\l'i,2uj'it) 2 

xs[u,w) = — — exp[{b+ -p{2uji3))wu], 

a[w/3\{2uJi/3, 2c<j3|)cr('u||2c<Ji/3, 2u!3\) 3 

xl{u,w/3)xl{u + 2uji/3,w/3)xl{u + 4uJi/3,w/3) , 

= 77^ 7^ 7^7 71 7^ 7^7 exp[bwu, 4.41 

xl[2uji/3,w/3)xl{4:Ui/3,w/3) 

VM = piu), Vsiu) = piu\{2uji/3,2uj^}). 

Here again, b is an arbitrary complex parameter. The short and long root functions also 
satisfy the A2 functional equation ( |4.32| ) separately, as expected. 

Finally, we consider the functional equation ( |4.26| ) for the l2{m) root system. The 2m 
roots are labelled in order of increasing angle by {ai, . . . , a2m} with A+ = {aj,j = 1, . . . , m} 
and aj^rn = for j = ^1 ■ ■ ■ For example, choosing all roots to have the same length 
\aj\ = 1, a possible basis is given in (2^). Unlike the crystallographic root systems, there is, 
in general, more than one functional equation. The functional equation (|4.26|) corresponding 
to = SpSa with ip = 2TTN/m, N = 1, . . . , [m/2] is 

d 



= /i ■ RN./m Q^fh(m) (4-42) 



with 



fh{m)i<l^ ^) = H 9\j\ 9\j+N\ x\j+N\{aj+N ■ q, aj+jv- • 9. -«J+2jV ^^)■ (4-43) 

In contrast to the previous cases, the the coupling constants are included in the functional 
equation since they do not factor out, in general. Because of the many functional equations 



16 



(|4.42| ) to be satisfied, only rational solutions are allowed in this case. For odd m, l2{fn) 



roots are all in a single orbit and only one coupling constant and function are possible. The 
solution to ([4.42|) in this case is 

xiu^w) = ( ")exp[6Mw]. (4.44) 

\u wj 

However, l2{fn) for even m has two orbits: one the set {aj} with odd j and the other with 
even j. It is possible to have two independent coupling constants and functions in this case. 
The corresponding functions are denoted xq and xe, respectively and the solution of (|4.42| ) 
is 

xo{u,w) = ( ^explbuw], (4.45) 

\u w J 

xe{u,w) = ( —^explbuw], 

\u wJ 

in which b and c are arbitrary complex constants. 

For the other non-crystallographic root systems H3 and H4, the two-dimensional sub-root 
systems are A2 and /2(5). Thus the solution ( |4.44| ) satisfies all of the functional equations 



for the consistency of the Lax pair. 

At the end of this section, let us show the Lax pair operator formulation for the rational 
potential with a confining harmonic force. This applies, as before, to all of the root systems 
including the non-crystallographic ones. This is a simple generalisation of the Lax pairs given 
in 0, which were constructed for the vector representations of the classical root systems, 

■Aj-^ Bj.^ Cf^ I-^T'i cLlld BCy. 

Let us start from the Lax operator for the rational potential without spectral parameter: 
L = p-H + iY. g\^\ {p'-H)x{p-q)Sp, (4.46) 

M = i 9\p\y{p-(l)Sp, 

1 1 
x{p-q) = , y{p-q) = -i ^ = -x{p-qf. (4.47) 

p-q [p-qY 

It corresponds to the following Hamiltonian and the equations of motion: 

Tr(L^) oc „.4p^+i:gi^. (4.48) 

L=|L,M1 « , = p = 2i:gi|^. (4.49) 
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If a confining harmonic potential is added to the Hamiltonian 



«.4/ + i^V+i:§^. (4.50) 



the equations of motion read 



,2 



g = p, p=-uj\ + 2Y.^^T^r (4-51) 



a|2 (a ■ g)3 



It is elementary to see that the above equations can be written in an operator form 

L=[L,M]-uj'^Q, Q=p-H = L-X, (4.52) 
in which L and M are the same as before ([4.46|) and Q is defined by 

Q = q-H. (4.53) 

It is easy to verify that 

[Q, M] = -X. (4.54) 

This property is valid only for the rational potential without spectral parameter. 
Next let us introduce a pair of non-hermitian operators L"^ by 

L^ = L±iLuQ. (4.55) 

Their time evolution equations read 

= L±iujQ 

= [L, M] - uj^Q ± iuj{L - X) 

= [L^ T iooQ, M] - u'^Q ± iu}{L - X) 

= [L^,M]-uj'^Q±iuj{L^ Ti^Q) 

= [L^,M]±iujL^. (4.56) 
If we define hermitian operators Ci and £2 

£i = L+L-, C2 = L-L+, (4.57) 

they satisfy Lax-like equations 

4 = [£„M], J = 1,2. (4.58) 
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From these we can construct conserved quantities 



Tr(4), j = l,2, A; = l,2,..., 



(4.59) 



as before. It is elementary to check that the first conserved quantities give the Hamiltonian 

Tr{Cj)<^n^, J = 1,2. (4.60) 

This then completes the derivation of the Lax pairs for all of the generalised Calogero-Moser 
models. 

5 Representations of the Lax operators 



We next consider representations of the operator algebra ( |4.1| - |4.4| ) for operators {Hj, j = 
l,...r} and {sa, ol G A} that appear in the Calogero-Moser Lax pair ( [4.5| ). In general, 
the representation of an algebra consists of a vector space V and a mapping from elements 
of the algebra to G'L(V), e.g. Sa lZ{sa). When a basis is chosen for V then lZ{sa) 
becomes a matrix in GL{d) where d is the dimension of V. First define a monomial Pn{t) = 



n 



,(i) 



,(J) . 



t) in an auxiliary vector variable t G associated with a set Q of vectors 



{v^^' G R^, J = 1, . . . A^}. The basis vectors of V, for a representation TZq of the algebra 



( |4.1H4.4| ), are monomials resulting from the orbit of Puit) under the reflection group, i.e. 
Pn{sa{t)) with a G AQ All of the monomials in V therefore have the same degree A^. A 
similar representation for only the reflection generators Sq with Q a sub-root system of A 
was introduced by MacDonald for irreducible representations of the Weyl groups [jll] , p!2 
The action of the operators on the basis vectors in V, is 



N 

i=i 

N 

'(V 



(i) . 



N 
N 



,0> 



t), 

N 



(5.1) 



Note that the representation matrices {TZ{Hi), / = 1, . . . , r} are all diagonal in this basis. 

All of the previous versions of Lax pairs for the Calogero-Moser models, namely the 
minimal and root type Lax pairs, result from ( ^.5| ) in a representation TZq with Q a single 

-'^More generally, the vector space for the representation is Sym{W*), the symmetric algebra on the 
r-dimensional dual vector space W* . 
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vector, i.e. (|5.1| ) with = 1 p]. In the case of = 1, the vectors in V, upon which the 
representation matrices act, may be denoted in a simple manner as with 

\fi)^{fi-t). (5.2) 

For the minimal type Lax pair, the representation has G V with /i being a weight in the 
same Weyl orbit as the highest weight of a minimal representation of the corresponding Lie 
algebra. For the root type Lax pair, the vectors in V are labelled by all of the roots for 
the simply-laced root systems and either the short or long roots for the non-simply laced 
root systems. All of the fi which label the basis vectors of V are then in a single Weyl orbit, 
as required. It will be demonstrated in section |] that the Lax pairs in these representations 
agree with those given in 0. 

There is a simple relation between the geometry of the root system and the dimension 
of the TZq representations for Q a single vector. The vectors fi labelling the \fi) G V are 
generated by the orbit of a single vector fiQ under the action of the finite reflection group. The 
vector fiQ may also be assumed, without loss of generality, to be in the principal Weyl chamber 
or its boundary. The dimension of the representation depends only on which reflection 
hyperp lanes, if any, contain the vector fiQ. Let $o be the set of indices j such that $o = 
{j, fiQ ■ aj = 0, aj G n}. For crystallographic root systems $o are the Dynkin labels of 
/io which are equal to zero. Then s^^. (/iq) = /io for j G $o or the isotropy group Jq of 
/io is generated by reflections about simple roots aj G 11 for which j G $o- Therefore 
the number of elements in the orbit of /iq under the reflection group W, and hence the 
dimension of the representation, is = |iy|/|/o|. Since the isotropy group is the direct 
product of the reflection groups corresponding to the connected parts of the Coxeter diagram 
for the original root system, after deleting those vertices corresponding to the indices j not 
contained in Jq, |/o| may be computed as the product of the orders of the reflection groups 
of the corresponding Coxeter sub-diagrams. 

As a first example, we calculate the dimensions of the minimal and root type repre- 
sentations of the algebra in ( [4.1| - [4.4[ ) and hence the Lax pairs for A^. There are three 
inequivalent minimal representations with /iq = A^^\ A^'^\ or A^^\ For /io = A^^-*, the vec- 
tor representation, the isotropy group Jq = Wa4, so D = \WAr^\/\WA4,\ = 6!/5! = 6. For 
the antisymmetric product of two vector representations, /io = A^^\ Iq = Wa^ x Wa^ so 
we obtain D = 6!/(2!4!) = 15. The antisymmetric product of three vector representa- 
tions, /io = A(^), Jo = Wa2 X Wa^ so we obtain D = 6!/(3!3!) = 40. Finally, for the root 
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representation, choose /iq = A*^^) + A^^\ which is the highest root. Then Iq = Wa^ and 
D = I 1 / 1 W^As I = 6!/4! = 30. As an aside, we note that the highest dimensional repre- 
sentation of the Lax pair for has /io = aiA^^'^ + a2A'^^^ + a^A^^^ + a^A^^^ + a^A^^^ with 
all non-vanishing coefficients aj,j = 1, . . . , 5. Since the dimension of the representation de- 
pends only on which Dynkin labels are non-zero, these coefficients aj may be chosen to be 
1, without loss of generality. Then the isotropy group is trivial and we obtain D = 6\, which 
has not been previously described. In fact, even among the representations with = 1, 
most of the matrices for the Lax pairs had not been described before. 

As a second example, let us evaluate the dimensions of some lower- dimensional Lax pairs 
for the Es model. Let us take fio = A^^^ and fiQ = A^''^ corresponding to the two end points of 
the two longer forks of the Eg Dynkin diagram. For /xq = A*^^^ the isotropy group Iq = We-j 
and D = | W^^s I / 1 I = 240. This gives the root type Lax pair. For /io = A^"^) the isotropy 
group Jo = Wd-j and D = \WEa\/\WD-j\ = 2160. This gives the second smallest Lax pair for 
the Eg model and its weights are a part of the 3875 representation. 

For the non-crystallographic root systems, the root type Lax pairs with dimensions 2m 
for him), 30 for and 120 for H4 give the smallest dimensional Lax pair matrices. 

6 Minimal and root type Lax pairs 

As examples of representations of the Lax operators, we consider those representations that 
yield the minimal and root type Lax pairs previously reported in Ref. 0. The functional 
equations associated with these representations are also derived by restricting those given 
in section ^ These are the same as those given in our previous papers, which we denote as 
sum rules. In this section we consider only the crystallographic root systems. 

Minimal type Lax pairs are associated with minimal representations of the Lie algebras. 
All of the fundamental representations of Ar, the spinor representation of Bj., the vector 
representation of Cr, and the vector, spinor and anti-spinor representations of Dj., the 27 
and 27 of Eg and the 56 of Ej are the minimal representations. All weights in a minimal 
representation are in a single Weyl orbit. The vectors G V of the representation are 
/U G in which T^mm is the set of weights of the minimal representation. They 

are characterised by the condition 

p^- /i = {0, ±1}, peA, HEJ^min- (6.1) 



21 



On the other hand the vectors G V of the root type Lax pairs are a G A} 

for a simply-laced root system A. The set of basis vectors of the two possible root type 
representations for a Lax pair associated with a non-simply laced root system are a G 
As}, the set of short roots, and {\a), a G A^}, the set of long roots. Let us collectively 
denote by A^ the set of basis vectors of various root type representations: 

{A, for simply laced models, all roots, 
Al, for non-simply laced models, long roots. 
As, for non-simply laced models, short roots. 

The roots are characterised by the condition 

p'^-a = {0, ±1, ±2}, p,aeA, (6.2) 

except for the G2 case in which ±3 are also possible. 

The fact that the eigenvalues of the operator p^- H are restricted to these values ( |6.1| ), 
(|6.2| ) simplifies the representation of the Lax pair operator ( |4.5| ) drastically. Especially at 
the zero eigenvalue of p^- H the functions x\p\ and y\p\ take the following simple forms: 

{—1/,^, untwisted, 
—2/^, twisted short roots, (6.3) 
—3/,^, twisted short roots, G2, 

lim y\p\ (u, w) = - V|p| (u) + DipI, (6.4) 

in which D\p\ are constants possibly depending on the orbits, that is, and Ds for the 
non-simply laced cases. 

With the basis vectors of V, upon which the representation matrices act, labelled by a 
single vector, as in (|5.2| ), the matrix elements of the minimal and root type Lax pairs are 

L\u) = Yl T^{L)^,\p), M\v)= 7^(M)^,|p), p,z/GE„i„, (6.5) 

rmn 

m = 5] 7^(L)„^|a), M\(3)= 5] 7^(M)„^|«), a,/3GA^, (6.6) 

where L and M are the Lax operators defined in ([4.5|) . Hereafter we adopt simplified notation 
Lpu {La/3) and M^^, {Ma/3) for the matrix elements 7l{L)p^ {7l{L)a/3 ) and 7^(M)^^ {7l{M)a/3 
). For an arbitrary function f{u) we have the following matrix element for the minimal type 

f{p''-H)sp\u)= Y {/(-1)Wp + /(0)5mAm-o + /(1)^™}|/^), (6-7) 

rmn 
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and for the root type (except for the G2 case) 



fip""- H)Sp\P) = {/(-2)^a-/3,-2p + /(-l)^.-/3,-p + /(0)5.A-,o 



a). 



(6.8) 



The G2 case can be handled in a similar way. 

By combining (|6.3| ), (|6.4j) and ( |6.7|) , (|6.8| ) it is straightforward to derive the matrix 
representations of the minimal and root type Lax pairs. For the minimal Lax pair we obtain 



P ■ At (5^,^ + i 9\p\x\p\ (p ■ q, i)5^,^u,p + ^m^p.i., 

p6A 



(6.9) 



pv 



E 9\p\V\p\{p ■ q) 

pGA,pVp=l 



peA 



in which is a constant and -B^ contains the dynamical variables q. Both have no effect 
on the Lax equation and can thus be omitted. They are 

= '-{gLNr + gsN^), for untwisted, '-{gLN^ + 2gsN^), for twisted, 

Bn. = UgLNTDL + gsm'Ds) - '-Y.9\p\V\p\{p- q), 

peA 

in which (^™) is the number of long (short) roots a G A such that a ■ /i = for a 
given /i G S^jn. The integer N"^ (^™) is well-defined since all vectors labelling the basis 
elements in V for the representation are in a single Weyl orbit and hence any other vector, 
say V = w{p) for w an element of the Weyl group. This implies that a ■ z/ = if and only 
if w{a) • yU = and so is the same for any choice of /i G ^min- When the and 
terms are dropped, the above Lax pair ( |6.9|) , ( |6.10| ) has the same form as that given in our 
previous paper 0]. 

For the root type Lax pair we obtain 



Laf3 = p-aSa^/3 + iJ29\p\[x\p\{p-q,0^a-(3,p + '^X\p\{p-q,2^)Sa-(3,2p 

peA 



M, 



af3 



9\a\V\a\{a ■ q) + E 9\p\y\p\{p ■ q) ^^,0 

A9p=Q:— (t,(t€ A 

+^ E 9\p\ [y\p\{p ■ q, 0^a-i3,p + y\p\{p ■ g, 2^)6a-f3,2p 

peA 



+ A,(5„,/3,(6.11) 
(6.12) 



+ Br 5 a, 3 ■ 
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As in the minimal case Ar is a constant and Br contains the dynamical variables q. Both 
have no effect on the Lax equation and can be omitted. They are 

Ar = ^{gLNl + gsN's), for untwisted, ^{gLN^ + 2gsN's), for twisted, 
Br = ligLNlDL + gsN^sDs)~lY.9Mp-q), 

in which {Ng) is the number of long (short) roots a E A such that a ■ (3 = for a 
given P G A, {A^, As). When the Ar and Br terms are dropped and with the following 
identification: 

xiiu,^) = xsiu,^) = x{u,^) or x^^/^\u,^), 

xl{u,2^) = xs{u,2^) = Xd{u,^) or x^];^'^\u,^), (6.13) 



the above Lax pair ( |6.11|) , ( |6.13|) has the same form as that given in our previous paper 



The restrictions on the eigenvalues of the operator ( |6.1|) , ( |6.2|) , simplify the functional 
equations, too. Let us examine the A2 functional equation (|4.32|) by adopting the variables 
u, V, ^1, and ^2 defined in (|A.8| ). We impose a condition ^1 = ^ = —^2, so that the restriction 



of the minimal representation (|6.1|) is satisfied. The limit formulas (|6.3|) and ( |6.4|) yield the 
so called first sum rule of Ref. 

x{u, y{v, - yiu, xiv, + x{u + v, [V{v) - V{u)] = O. (6.14) 

Here the suffix |p| is omitted since all the roots participating in the A2 functional equation 
belong to the same orbit. Likewise, let us impose a condition ^1 = 2^, ^2 = so that the 
restriction of the root type representation { p.'2\ ) is satisfied. The limit formulas ( |6.3| ) and 
yield the so called second sum rule of Ref. P| 



= x{u + v,^)y{u,^)-y{u + v,^)x{u,^) 

+2 [xiu, 2^)y{v, - vi-v, Ox{u + v, 2^] 

+xi-v, Oviu + V, 20 - yiu, 2^)xiv, 0- (6.15) 

One may also verify that all of the functional equations in which must be satisfied 
by the functions x\p\{u,w) appear as various restrictions of the operator equation in (|4.27|) , 



although we will not derive them here. Since the functional equations for the minimal and 
root type representations are restricted, their solution spaces are larger than that for a 
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generic representation. For example, these functional equations have the same solutions as 
the general solutions given in section ^, except that the exponential factor is changed from 
exp[bwu] to exp[b{w) u], in which an arbitrary function b{w) need not be linear in w. 



7 Comments and Discussion 

Firstly, let us comment on the relation between our work and the paper on elliptic Dunkl op- 
erators by Buchstaber, Felder, and Veselov ^ . From the commutativity of Dunkl operators 
they derived a functional equation Eq. (9) of Ref. [Q, which was closely related to our equa- 
tions = 0, = 0, = 0. But their equation did not contain the spectral parameter 
dependence explicitly. They obtained what would amount to our untwisted solution for the 
/yij = functional equation. This gives a clue that the classical and quantum integrability 
of the generalised Calogero-Moser models are closely connected. Secondly, some remarks 
about the Calogero-Moser models based on Br, Cr, and BCr root systems. The short roots 
of Br, the long roots of Cr, and the long and short roots of BCr, when restricted to any 
two-dimensional plane, form only an Ai x Ai sub-root system. This means that their short 
root function xs{u, w) (and/or xl{u, w)) in these models are required to satisfy the B2 func- 
tional equation only but not the A2 one. Thus the solution space is larger than that of the 
other models. This in turn allows more potential function terms (one more for Br, and Cr 
and two more for BCr) with independent coupling constant (s) in the Hamiltonian without 
breaking integrability. We call these models extended twisted models. The root type and 
minimal type Lax pairs for the extended models are given in |T^ . The Lax pair operators 



for the extended twisted models can be constructed in a similar way as is given in this paper. 
Thirdly, a few words about the structure of the functions x\p\{p ■ q, {p^- /i)^) and their func- 
tional equations ( [4. 261 ), ( |4.27| ), and the self-duality of the two-dimensional crystallographic 



root systems. The coefficients in the functional equation ( [4. 261 ), i.e. p^- fi and Sp{aY- p, etc., 
come from the second argument of the function x\p\{p-q, (p^- p)C)- Namely, they are co-roots. 
On the other hand, the gradient operator, d/dq, in (|4.27| ) supplies the coefficient from the 
first argument of the function x\p\{p ■ q, (p^- p)0) that is the roots. The operator -R±,r/4 in 
the B2 case and R±Tt/6 in the G2 case rotates a short root into a long root position and vice 
versa. In other words, the rotation operator R^/2 in (f4.27 ) performs the necessary conversion 



from the roots to the co-roots. This is possible because of the well-known self-duality of the 
two-dimensional crystallographic root systems, A2, B2, and G2, under a = lajo?. 
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As a final remark, let us comment on the integrability of the generalised Calogero-Moser 
models. As is well known, the existence of the independent involutive conserved quantities 
as many as the degrees of freedom is necessary and sufficient for the Liouville integrability. 
To the best of our knowledge, the involution of the conserved quantities for the models based 
on the exceptional root systems and the non-crystallographic root systems as well as all the 
twisted models is yet to be demonstrated. 
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A Appendix: Solutions to functional equations 



In this appendix we demonstrate that the solutions (|0^ ), (|07| ), (|Og| ), P^ ), 



( |4.45| ) given in section ^ satisfy the functional equations for the consistency of the Lax pair 

d 

= —f{q,fi). (A.l) 

We do not give a proof that these solutions are the most general ones. First let us remark 
on a symmetry of the solutions of ( [A.l|) such that if x\p\{u,w) is a solution then 



x'^p^{u,w) = x\p\{u,w)exp[buw], V6 G C, (A. 2) 



is also a solution. Therefore, it is necessary to prove ( [A.7| ) only for one representative function 
x\p\{u,w) among those related by the symmetry (|A.2| ). Assume that the function x\p\{u,w) 
satisfies (|A.1|) with f{q,p) defined by ([4.28|) . The corresponding expression in ( [4.28|) for the 



new solution xl^p<^^{u,w) = x\p\{u,w) exp[buw] will be denoted as /(g, /i)"^"". Then 



X exp [b^ {{p . q){p^. /i) + (a ■ p))] \. (A.3) 
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The exponent is proportional to 

q- fx-q- (Spifi)) = q- fx-q- R-^(X. 



Two pairs of roots (p, a) and (p', a') which appear in different terms of ( |A.3D are related by 
SpSa- = Spis^i, which implies also 

s^Sp = s^iSpi = R_^. (A.4) 

Therefore the exponential factor in ( |A.3| ) is common to all terms and may be factored out 
of the sum. The functional equation ( [A. ID for /((?, p)"'^"' then is 



^^■R^,2 ^-^fiq^f^r"" (A.5) 
d 

= ■ ^ (/(g, p) exp [b^{q ■ fi-q- R^^^i)]) 

= b^fiQ, At) exp [b^{q ■ p - g ■ R-^pjx)] (p • R^/2l^ - At ■ i?-v,/2/i) 
= 0. 

Here the orthogonality of rotation operators 

P ■ -Rv/sAt ~ At ■ R-ip/2fJ' = 0, R-ip/2 = (^-6) 

is used. This then demonstrates that the transformed function x\p\{u,w)^'^^ satisfies the 
functional equation if the original function x\p\{u,w) does. 

Next we show in turn that all of the solutions given in section ^ satisfy 

/G(g,p)=0, G = A2,B2,G2, (A.7) 

which is a sufficient condition for the solutions of ( A.l ). Note that this is a necessary and 



sufficient condition for the commutativity of the corresponding Dunkl operators 
A.l A2 solution (i^) 



Defining the variables 

u = a ■ q, V = P ■ q, (A. 8) 

^1 = a^-p, ^2 = /?^-/U, 

/yi2(g, p) becomes 

/A2(g,/^) = X{u + V,^i+ ^2)x{u,-^2) +X{v,^2)x{u + V,^i) - x{u, ^i) x{v, + ^2)- (A.9) 
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By substituting the solution 



a(w — u) / ■ X 

x{u,w)= ; ; A.IO 



it reads, after extracting a common denominator, 

fA,{u,v,^u^2) = [(r{^i + ^2-u-v)a{^2 + u)a{v)a{^i) (A. 11) 

- v) cr(6 -u-v) a{u) a{^i + ^2) 

- u) a{^i + ^2-v) aiu + v) ^(6)] 

/ [aiu) aiv) aiu + v) ^(^1)^(6) ^(6 + 6)] • 
The expression in the numerator may be shown to vanish, by using the following identity 



(see page 153 of and Eq. (5.1) of ^ 

= (7(5; - ui) o-(2; + ui) (j(u2 - U3) (j(u2 + M3) (A.12) 
+aiz - U2) aiz + U2) aius - Ui) aius + Ui) 
+aiz - U3) aiz + U3) cr(Mi - U2) cr(Mi + U2) 

with the choice of variables 

z = 6 + ^(6-^), ui=u-^i^,-v), (A.13) 

U2 = ^i^l + v), U3=^i^i-V). 

This implies that fA_2 iQ, /^) = for the function in ( |A.1CI|) and after including the symmetry 
of solutions ( |A.2| ), the more general function in ( [1.34| ) also gives /^al?''/^) = therefore 
the A2 functional equation is satisfied. 

A. 2 B2 untwisted solution (|4.csV|) 

Next we turn to the functional equation arising from the B2 sub-root system. This, as well 
as the functional equation associated with the G2 root system, admits two types of solutions, 
the untwisted and the twisted ones. By using the same definitions of the variables u, v, ^1, 
and ^2 as in (|A.8|) , except for the fact that a and f3 are the short and long simple roots of 
B2, respectively, the expression for fB^ililA becomes 

fB2iu,v,ii,i2) = - xsiu,ii)xLiv,ii+ ^2) (A.14) 
+ xsiu + v,^i + 2i2)xLi2u + v,-i2) 
+ xsiu, -^1 - 2^2) XLi2u + v,ii + 6) 

+ Xsiu + V,^i)xLiv,i2)- 
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Substituting in the untwisted solutions 



xl{u, w) = xs{u, w) = ^\ (A.15) 

(r[w)(J[u) 



this becomes 



fB,{u,v,^i,^2) = [-(r{^i-u)a{^i + ^2-v)a{u + v)a{2u + v)a{^2)cr{^i + 2^2) 

+ 2^2-u-v) a(6 + 2u + v) a{u) a{v) a{^,) (1(^1 + 6) 
+^(6 + 26 + u) (t(6 + ^2-2u-v) aiv) a{u + v) a(6) a(6) 
+^(6 - u - t;) a(6 - t;) (t(2u + v) a(6 + 6) ^(6 + 26)] 
/ [o-(^i) a{v)a{2u + v) a(6) a(6) ^(6 + 6) ^(6 + 26)] • (A.16) 

Let us denote the numerator by (7b2(^) ^? 6? 6)- Gathering terms, one obtains 

5'B2(w,w,6,6) = or(2M + t;)cx(6 + 26)[-o-(6 -m)o-(6 + 6 -^^)o-(w + ^^)cr(6) 

+a(6 - u - t;) a(6 - v) aiu) (t(6 + 6)] 

+(r{v) fT(6)[a(6 + 26 - M - 1^) o-(6 + 2u + v) a{u) a(6 + 6) 

+(x(6 + 26 + m) ^(6 +^2-2u-v) a{u + t;) (t(6)]. (A.17) 



Using the identity (|A.12|) for the expression in the first set of brackets with the choice of 
variables 

z = -i^ + ^(-M + 6 + 6), «i = ^(n-6 + 6), (A.I8) 
U2 = ^(^^ + 6 + 6), Ms = ^(-u + 6 + 6), 

and for the expression in the second set of brackets with the variables 

z = ^(m + 6 + 36), mi = ^(3m + 2t;-6-6), (A.19) 

U2 = ^(m + 6 + 6), m3 = ^(-m + 6 + 6), 



gives 



5(^2(^,^,6,6) = -o-(2u + t;) a(6 + 26) o-(?;) a(6 + 6 - « - ^) o-(6) cr(6 + m) 
+a{v) a(6) a(6 + u) (t(6 + 26) a(6 + ^2 ~ u - v) a{2u + v) 
= 0. (A.20) 
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Including the possible symmetry transformations ( |A.2| ) of the function (|A.15|) this then 



proves that /bjIq')/^) = for the general elliptic untwisted solution of (|4.37|) . 

The untwisted solution ( [4.37|) of the G2 functional equation ( [4 .391) may be proven again 
using only the a function identity ( [A.12|) . Since the method of proof is essentially the same 



as for the B2 functional equation we omit the details of the proof. 
A. 3 B2 twisted solution (i^) 

We next demonstrate that the twisted solution in ([4 .381) satisfies the B2 functional equation 
( |4.35| ). First we define a particular untwisted solution Xq{u, w) of ( |4.37|) by assigning a special 



value of the constant b = rji/ui 

xs{u,w) = xl{u,w) = xo{u,w) = ^1^, exp[{r]i/uji)wu], (A.21) 

a[w)a{u) 

where r]i is defined in terms of the Weierstrass ( function as rji = C(^i)- The value of b is 
chosen so that xo{u,w) is periodic in the ui direction 

Xo{u + 2u!i,w) = Xo{u,w). (A. 22) 

Adding a constant vector to q does not affect the equation fs^ {q, /x) = so g is shifted as 

q^q+ , A.23 

in which A*^") is the fundamental weight dual to the simple short root a, i.e. 

a^-A(") = l, /?^-A(")=0. (A.24) 

Then the equation fs^ {q, /^) = reads 

0= — Xo{a ■ q + uija"^ ■ fi) xo{P ■ q, {2a + py ■ ^) 

+ Xo{{a + /3) ■ q + LUi,{a + (3y ■ fi) Xo{{2a + (3) ■ q, -0^ ■ /i) 

+ Xo(a ■ q + uJi, -{a + I3Y ■ /i) Xo((2q; + (3) ■ q, {2a + (3y ■ fi) 

+ Xo{{a + (3)- q + uji,a^ ■ ^)xo{(3 ■ q,(3^ ■ n), (A.25) 

where the periodicity of xq{u,w), ( |A.22| ), is used. This simply means that 

a{w-u) 
xl{u,w) = exp[bwu\, 
a[w)a[u) 

a{w — u — uJi) 

xs[u,w) = —. — ^ — 7 r expmi w + oiy M , (A. 26) 

^ ' o{w)o{u^ujx) 
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satisfy the equation /bjIq')/^) = 0. Since xl{u,w) here is the same function as in (|4.37|) 



any hnear combinations of xs{u, w) from the untwisted solution in ( [4.37 ) and in ( A.26|) also 



satisfy the B2 functional equation. Requiring that the linear combination should satisfy the 
A2 functional equation and that it has a simple pole with unit residue at m = 0, we obtain 



aiw — uj aiw — u — uJij 

exp[riiw\ 



exp[bwu]. (A. 27) 



This expression for x^{u,w) has the following monodromies in u and w: 

x^g{u + uji,w) = x^g{u,w) exp[(—r]i + buji)w] , 

x^g{u + 2lj3,w) = Xg{u,w) exp[2{—ri3 + b 003)10] , 

x^g{u,w + 2ui) = Xg{u,w) exp[2{—rii + bui)u] , 

w + 4CJ3) = Xs{u,w)exp[A{-r]3 + buj3)u]. (A.28) 

It also has the following poles and zeros in the fundamental regions of u and w: simple poles 
at M = and w = with residues 1 and —2, respectively, and a zero at u = w/2. It may 
also be shown that the twisted solution to the B2 functional equation xs{u, w) in ( [4.38| ) has 
the same monodromies and poles. This implies that the ratio x^g{u, w)/xs{u, w) is an elliptic 
function in both u and w and has no poles and therefore is a constant. Since the residues 
are the same at all the poles, the ratio is equal to one and x'g{u,w) is, in fact, the same as 
the twisted solution in ( |4.38| ) 

A.4 G2 twisted solution ( |4.41| ) 



The proof that the twisted solutions ( |4.41|) to the G2 functional equation ( |4.39|) follow from 



the untwisted solutions proceeds in a similar manner. We start from the particular untwisted 
solution xs{u, w) = xl{u, w) = xo{u, w) satisfying /g2(q'5 /^) = 0, with the periodicity ( |A.22|) , 
and shift q as 

q^q + ^^^, A.29) 

in which A*^"") is the fundamental weight dual to the simple short root a. We obtain 

= Xo{a ■ q+'^,a'^ ■ n)xo{(3 ■ q,{3a + (3^- H) 

+Xo{a ■ q + -(2a + /?)^- fi) xo((3a + /?) ■ g, (3a + py ■ /i) 
+a;o((2a + /3) ■ g + (2a + py ■ /i) a;o((3a + (3) ■ q, -(3a + 2/3^ ■ fi) 
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+xo((2a + P)-q + ^,-{a + {Sf ■ ^i) Xo((3a + 2/3) ■ g, (3a + 2/5)^- /i) 
+a;o((a + /?) ■ g + (a + /i) Xo((3« + 2/?) ■ g, /i) 
+a;o((a + /3) ■ g + a^- yu) xo(/3 ■ g, Z?"^- /x). 



(A.30) 



Shifting g again by the same amount, or in one step 



g ^ g + 



3|a| 



(A.31) 



gives the following equation 
= 



xo(a • g H — r^, a^- yu) a;o(/5 ■ 9; (3a + ■ At) 



3 

+a;o(Q! ■ g + -(2a + PY ■ fi) Xo((3a + /?) ■ g, (3a + ■ ix) 
+Xo((2a + /?) ■ g + (2a + /3)^- yu) Xo((3a + /5) ■ g, -(3a + 2/5)^- /x) 
+xo((2a + /?) ■ g + -(a + /5)''- /i) xo((3a + 2/3) ■ g, (3a + 2/3)^- yu) 
+Xo((a + /?) • g + (a + /5)^- /i) xo((3a + 2/?) ■ g, -/5^- /i) 
+a;o((a + (3) ■ q + a"^- /i) Xo(/5 -q^P^- n), 



(A.32) 



where the periodicity of xo{u, w), ( \A^ , is also used. Adding ( |OoD and ( [02D together and 
using the symmetry of solutions in ( |A.2| ) implies that the following functions are solutions 
to the G2 functional equation (|4.39| ): 



xl{u,w) 



xs{u,w) 



a(w — u] 



a{w)a{u) 



exp[bwu] 



aiw)aiu+^-^) ^^P[^2/3),, w] + ^^^^^^^ _^ 4^, exp[(4/3)r,, w] 



(A.33) 
exp[bwu]. 



3 J a{w)a{u + 

Again, since xl{u, w) is the same function as in ( [4.37| ) an arbitrary linear combination of the 
xs{u,w) from ( |4.37| ) and ( [A.33| ) also satisfies the G2 functional equation. Requiring that the 
linear combinations should satisfy the A2 functional equation and that it has a simple pole 
with unit residue at m = 0, we obtain 



Xs{u,w) 



aiw — u) 



+ 



a[w — u — 



2uJi ■ 



— exp[(2/3)r7i w] 



{AM) 



a{w)a{u + 



exp[(4/3)r7i w] 



exp \pwu\. 
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This function has the following monodromies: 



x^g{u, w) exp 



-r]i + buji)w 
x°s{u,w) exp [2(-?73 + bu3)w] 
x° (m, w) exp [2(— ?7i + bui)u] , 
Xs{u,w) exp [6(-?73 + bu3)u] . 



(A.35) 



The following properties of x^(-u, w) can be shown. It has poles and zeros in the fundamental 
regions of u and w. simple poles at m = and w = with residues 1 and —3, respectively, 
and a zero at m = w/3. Since it may be shown that the twisted solution xs{u,w) to the G2 
functional equation in (|4.41| ) has the same monodromies, poles and zeros as x^g{u, w), by the 
same argument as in the B2 case given above, x'g{u, w) = xs{u, w) and the solutions derived 
above are the same as the twisted solutions in (|4.41| ). 

A. 5 Dihedral solutions (i^), (i^) 

Next, it will be shown that the solutions ( [4.44| ) and ( [4 .451) satisfy the dihedral l2{iTi) func- 



tional equations (442) for odd or even m, respectively. We assume that all roots are of 
the same length, even for an even integer m, since they may be made so by a redefini- 
tion of the coupling constants. First consider the case of an odd integer m and arbitrary 
N = 1, ... ,[171/2]. Substituting the functions 

x(u, w) = ( ^ 

\u wj 

in ([4.43|) for //^(m) (q', yu) and redefining — > ^|ajp/2 gives 



(A.36) 



pN 

J him) 



+ 



1 1 



+ 



{aj+N ■ q) {aj ■ q) {aj+N ■ q) (aj+2jv ■ /^)^ 
1 1 



The 0(1/^) terms cancel pairwise using the identity 

j 
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E 



,■=1 (ttj+d ■ q)iaj+N+d ■ Ai) ,=1 (ttj 



+d+h 



q){aj 



j+d+N+h 



(A.37) 



(A.38) 



with h, d, and arbitrary integers. This identity follows from the property of the l2{'m) 
—cxj- The vanishing of the and 0(1/^^) terms follows from a simple 



roots a 



trigonometric identity (g = |g| (cost, sin i:)): 



0. 



(A.39) 



~[ cos(t — jir/m) cos(t — (j + N)^ /m) 

The left hand side is a meromorphic function in t with a period vr and it is exponentially 
decreasing as t ^ ±ioo. It is elementary to show that all the residues of the possible simple 
poles f + J = 1,2, ... ,m vanish. Thus ( |A.39D vanishes. This shows that the dihedral 
functional equation ( |4.42| ) for odd m is satisfied by the functions in ( |A.36| ). Including the 
symmetry of solutions of ( A. 2 ) implies that the more general solution in (^]4^) also satisfies 



(ill). 



Next consider the case of even m and odd A^ in the functional equation ([4.42|) . Substi- 
tuting the functions 

xo{u, ly) = f- - — y xe{u, w)= {-- — ) 
\u wj \u wj 

in the equation for /^(^^(q', /i) ( [4.43D and redefining ^ ^\aj\'^/2, as before, gives 



(A.40) 



m/2 

fh{m)iQ^ /^) = 9o9e E ^E{(y2j ■ q, -(a2i+27v ■ xoia2j+N ■ q, ia2j+N ■ 

i=i 

+XE{a2j+N-i ■ g, {a2j+N-i ■ xo{a2j~i ■ q, -{a2j+2N-i ■ 1^)0 
1 c 



9o9e\ J2 

{j=l 



_{aj ■ q){aj+N ■ q) {.(^j+2N ■ fJ'){oij+N ■ At)^^ 

m/2 



1 f^i[ioi2j+N ■ q){o:2j+2N ■ IJ) (tt2i-i ■ q){a2j+N-\ ■ A*) 



+ 



{a2j ■ q) {a2j+N ■ At) («2j+Af-i ■ q) ("2, 



\j+2N-l ■ 



/i) 



(A.41) 



As before, the 0(1/^) terms in fjl(^rn)iq^ f^) cancel pairwise using the identity 



1 



m/2 , m/2 

E = E 

j^l {a2j+d ■ q){oi2j+N+d ■ /^) j=i {0l2j+d+h ■ q){oi2j+d+N+h ' /^) 



(A.42) 



with h an even integer and d and A^ arbitrary integers. The 0(^°) and 0(1/^^) terms 
are proportional to the corresponding terms for m odd given above and therefore vanish. 
Including the symmetry of ( |A.2[ ), the solution of ( |4.45| ) satisfies ( [4.42| ) for m even and A^ 
odd. 



34 



For the case of the solutions for the l2{'m) functional equation for even m and even iV, 
note that {q, yu) niay be written as 

m/2 



goXo{a2j+N ■ q, {a2j+N ■ 1^)0 xo{a2j ■ g, -{a2j+2N ■ 
+gEXE{a2j+N-i ■ <?, {a2j+N-i ■ a;£;(a2j-i ■ q, -{a2j+2N-i ■ 



(A.43) 



Here gfjijn/2)i^^ f^) proportional to /^) with the /2(m/2) roots rotated by vr/ 

and so the vanishing of //^(m) (O'; /^) follows using the previous equation and induction on A^. 
Therefore the solutions in ( [4 .451 ) solve the l2{iTi) functional equation ( |4.42| ) for all m and A^. 



m 
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